Mon. Not. R. Astron. Soc. OOO.[THT7lC2010') 



Printed 12 July 2011 (MN WF^ style file v2.2) 



X-ray flare candidates in short gamma-ray bursts 

R. Margutti^'^ *, G. Chincarini^'^, J. Granot^''^'^, C. Guidorzi^, E. Berger^ 



O 



oo 



Oh- 

o 



> 

o 



o 



X 



M.G. Bernardini^, N. Gehrels 

^ INAF Osservatorio Astronomico di Brera, 
^ Univerisitd Milano Bicocca, Dip. Fisica G. 



, A.M. Soderberg'', 

via Bianchi 46, Merate 23807, Italy 
Occhialini, P.zza della Scienza 3, Milano 20126: 



M. Stamatikos^, E. Zaninoni^'^^ 



Italy 



Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 

Raymond and Beverly Sackler School of Physics & Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 
Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hateld, Herts, AL109AB, UK 
University of Ferrara, Physics Dept., via Saragat 1, 1-44^'^'^ Ferarra, Italy 
Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA02138 
NASA-Goddard Space Flight Center, Greenbelt, Maryland 20771 

Dept. of Physics and Center for Cosmology and Astro-Particle Physics, Ohio State University, Columbus, OH 43210, USA 
' Univerisitd di Padova, Dip. Astronomia, v. dell' Osservatorio 31, Padova 35122, Italy 



Accepted Year Month Day. Received Year Month Day; in original form Year Month Day 



ABSTRACT 

We present the first systematic study of X-ray fiare candidates in short gamma-ray 
bursts (SGRBs) exploiting the large 6-year Swift database with the aim to constrain 
the physical nature of such fluctuations. We find that flare candidates appear in differ- 
ent types of SGRB host galaxy environments and show no clear correlation with the 
X-ray afterglow lifetime; flare candidates are detected both in SGRBs with a bright 
extended emission in the soft 7-rays and in SGRBs which do not show such compo- 
nent. We furthermore show that SGRB X-ray flare candidates only partially share the 
set of observational properties of long GRB (LGRB) flares. In particular, the main 
parameter driving the duration evolution of X-ray variability episodes in both classes 
is found to be the elapsed time from the explosion, with very limited dependence on 
the different progenitors, environments, central engine life-times, prompt variability 
time-scales and energy budgets. On the contrary, SGRB flare candidates signiflcantly 
differ from LGRB flares in terms of peak luminosity, isotropic energy, flare-to-prompt 
luminosity ratio and relative variability flux. However, these differences disappear 
when the central engine time-scales and energy budget are accounted for, suggesting 
that (i) flare candidates and prompt pulses in SGRBs likely have a common origin; 
(ii) similar dissipation and/or emission mechanisms are responsible for the prompt 
and flare emission in long and short GRBs, with SGRBs being less energetic albeit 
faster evolving versions of the long class. Finally, we show that in strict analogy to the 
SGRB prompt emission, flares candidates fall off the lag-luminosity relation defined 
by LGRBs, thus strengthening the SGRB flare-prompt pulse connection. 

Key words: gamma-ray: bursts - radiation mechanism: non-thermal -X-rays 



1 INTRODUCTION 

With an isotropic peak luminosity up to lO^^ergs"^, 
gamma-ray bursts (GRBs) are the brightest objects in the 
7-ray sky during their short lives (At ~ 0.1 — 100 s). 
Their duration-spectral hardness distribution gives evidenc e 
for the presence of two classes (|Kouveliotou et al.l 119931 ) : 
long and short GRBs (observed duration longer and shorter 
than 2 s, respectively), with short bursts appearing slightly 
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harder. The dichotomy in the duration-hardness dimensions 
suggested separate progenitor populations. However, until 
a few years ago, the distances, energy and environments of 
SGRBs (short GRBs) remained highly uncertain due to the 
poor localisation. 

The breakthrough in the study of SGRBs occurred 
thanks to the rapid slew capabilities of the Swift spacecraft 
(jGehrels et al.l 120041 ) which allowed spectroscopic observa- 
tions to be performed at very early times. These observa- 
tions revealed that SGRBs are cosmological, with prompt 
luminosities comparable to LGRBs albeit significantly less 
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energetic; with similar afterglows (jNysewander et al.ll2009l ) 
but residing in completely different environments. In sharp 
contrast to LGRBs, short bursts have been localised both in 
early-type and late-type host galaxies (see lBergeilbOlfl and 
references therein), pointing to an old progenitor popula- 
tio n. The detection of supernovae associated to LGRBs (see 



e^ 



„ iKulkarniet al.|[l993 : IStanek et aLll2003l : iFruchter et all 
l200a and references therein) provided instead support to 
models invoking young stellar progenitors. According to 
the standard scenario, LGRBs orig inate from the collapse 
of ra pidly-rotating, massive stars (jMacFadven fc WooslevI 
1 1999 ). while SGRBs are believed to result from the coales- 
cence of a binary system of compact objects (neutron star 
plus neu tron star NS-I-NS or neutron star plus black hole 
NS+BH, |Paczvnski|[T98^ : lEichler et alll 19891 : iNaravan etHI 
1 19921 ). 

Despite fundamental theoretical and observational 
progress, the nature of SGRB progenitors remains elu- 
sive. Numerical simulations show that the active stage of 
a NS-I-NS merger typically lasts ^ (0.01 — 0.1) s (see e.g. 
iNakar et al.ll2007l and references thereiqj): material ejected 
during the merger is expected to accrete on time-scales of 
the order of 1 — 10 s (the exact value depending on the ac- 
creting disk viscosity parameter and details of the ejection 
process). Thus, the detection of central engine activity on 
time-scales much longer than the usual dynamical or even 
viscous time -scales would chal lenge the currently accepted 
scenario (see lNakar et al.ll2007l for a recent review). 

Long-lasting (At ^ 10 s), soft energy tails detected in 
several SGRBs during t heir prompt emissio n (the so-called 
extended emission, see iNorris et al.l l2010al and references 
therein) represent such a case and pose severe constraints 
to existing models, especially when ene rgetically dom inat- 
ing with respect to the primary burst (jPerlevI |2009| V The 
same is true for the rec ently discovered presence of precur- 
sors l|Troia et al.ll2010h . Equally challenging would be the 
detection of late-time central engine activity in the form of 
flares superimposed over the smooth SGRB X-ray afterglow. 

Flares are currently detected in ^ 30% of long 
GRBs X-ray afterglows (jChincari ni et al.l [2010|) as fast- 
rise exponential-decay features whose spectral and temporal 
properties have been dem onstrated to show a str ict analogy 
to LGRB prompt pulses (JMargutti et al.ll2010bh : this find- 
ing suggested that flares might originate from re-activations 
of the LGRB central engine. Several ideas on how to explain 
the possible presence of flares in short GRBs have been ex- 
plored as well: the fragmentation of the outer parts of an 
hyper-accreting disk around the newly formed black hole as 
a result of gravitational instabilities could potentially lead 
to large-amplitude variations of the central eng ine output 
of both long and short GRBs jPerna et aLlbOod l. Alterna- 
tively, the late-time accretion of material launched into ec- 
centric but gravitationally bound orbits during the compact 
binary merger could provide the fuel to revive the central 
engine activity (IRosswog 2007l l. The long term evolution of 
debris following t he tidal disruptio n of compact objects has 
been identified bv lLee et al.l l|2009f ) as a feasible mechanism 
to produce fiares. Finally, as an alternative in the context of 
accretion-powered models, magnetic halting may also give 



^ In a recent study iRezzoIIa et al.l 1 12011^ found Af ~ 0.3 s. 



ris e to secondary episode s of delayed activity as suggested 
bv lProga fc Zhand i|2006r ). However, the observational prop- 
erties of fiares in SGRBs have not been determined, yet, so 
that it is at the moment unclear if any of these models would 
be able to explain the observations. 

While SGRB X-ray light-curves clearly show temporal 
variability superimposed over a smooth decay, the presence 
of real fiares in short bursts is questionable. In particular, it 
is at the moment unclear if what is currently identified as 
short GRB fiare emission (see e.g. iLa Parola et al.ll2006l for 
GRB 051210) quantitatively shares the very same properties 
of the population of long GRB flares: are there fast varying 
Af/f <ti 1, prominent temporal features in the afterglow of 
SGRBs with properties reminiscent of the long GRB flaring 
emission? Do SGRB flare candidates follow the entire set of 
relations found from the analysis of real fiares in long bursts? 
In particular: is the evolution of their temporal and energetic 
pr operties compatib l e with the flare-like behaviour identified 
by IChincarini et al.l (|201Q )? What is the typical amount of 
energy released during such episodes of variability? Is there 
any link between the late-time variability which appears in 
the X-ray afterglow of SGRBs and their prompt emission? 
Negligible spectral lag is a defining characteristic of SGRB 
prompt pulses: is this picture still valid when considering 
their late-time variability? 

Prompted by this set of still open questions, we present 
the first systematic study of X-ray flare candidates in short 
GRBs, taking advantage from the large Swifl 6-year data- 
base. Through a homogeneous temporal and spectral anal- 
ysis of the widest sample of SGRB light-curves available at 
the time of writing, this study allows us to perform a one-to- 
one compariso n with the properties of X-ray fiares detected 
in long bursts JChincarini et al.ll2010l.lMargutti et ahllJOlObl . 



iBernardini et al.ll201ll . iMargutti et al.ll201ll '). The primary 
goal of this paper is to observationally constrain the origin 
of SGRB flare candidates providing the reader with a com- 
plete picture of their properties. 

This work is organised as follows: the sample selection 
and data reduction is presented in Section[2] Results are de- 
scribed in Section [3] and discussed in Section 2] Conclusions 
are drawn in Section [5] 

The GRB phenomenology is presented in the observer 
frame. Isotropic equivalent luminosities and energies are 
listed. The observer frame 0.3-10 keV energy band is adopted 
unless specified. The zero time is assumed to be the trig- 
ger time. We use the notation: y|C!RB |-yLGRB-| ^^ j^j-_ 

cate that Y refers to the fiare (prompt) emission of SGRBs 
(LGRBs) . All the quoted uncertainties are given at 68% con- 
fidence level (c.L). Standard cosmological quantities have 
been adopted: Ho = 70Kms"^ Mpc"\ Oa = 0.7, ^m = 0.3. 



2 SAMPLE SELECTION AND DATA 
ANALYSIS 

We select the sho rt GRBs detected by th e Swifl Burst Alert 
Telescope (BAT. iBarthelmy et al.ll2005al ) an d promptly re- 
point ed by the Swifl, X-Ray Telescope (XRT. lBurrows et al.l 
I2OO5I ) between April 2005 and February 2011. The short na- 
ture of each event is established using the combined infor- 
mation from the duration, hardness and spectral lag of its 
prompt 7-ray emission: a prompt 7-ray duration Tgo ^ 2 s 
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coupled to a hard 7-ray emission with photon index F < 1.5 
and a negligible 7-ray spectral lag t^ are considered in- 
dicative of a short GRB nature (see Table [l]). The mor- 
phology of the host galaxy is also used as an additional 
indicator, when available. The final sample comprises 60 
SGRBs. The presence of X-ray vari ability in each SGRB i s 
investigated following the method bv lMargutti et al.l l|2011f ). 
used to determine the presence of flares in long GRBs. Only 
GRBs showing fluctuations with a minimum 2o|j signif- 
icance with respect to the continuum have been consid- 
ered in the following analysis. This procedure automati- 
cally identifies the best time intervals to be searched for 
the presence of X-ray flare candidates in SGRBs. Out 
of ~60 Swift SGRBs, 8 satisfy the variability require- 
ment above (Table [T]jj. Notably, the sample includes the 
unique 2 SGRBs with secure early-type host identification: 
GRB 050724 (Barthclmy et al. 2005b) and GRB 100117A 
iFong et al.ll20 10). In three cases (GRB 050724, GRB 070724 
and GRB 071227, in boldface in Table[T]) an extended emis- 
sion (EE) has been detected in the soft gamma-ray en- 
ergy rang e after the short hard spike iJNorris et al.ll2010al : 
iNorris et al. 2011). In the other cases, an upper limit on the 
EE to IPC (Initial Pulse Complex) i ntensity ratio (f?int = 
EEint/IPCint) has been provided by iNorris et al.l (|2010al '): 
for the sample of events without EE the upper limit on 
Rint is found to be a factor > 10 below the typical 7?int 
of SGRBs with detected EE (Table [T] column 7). Finally, 
GRB 100816A has not been included in the s ample in spite 
of its T90 = 2.9 ±0.6 s iMarkwardt et al.llioTol ') since the low 
statistics prevents the 7-ray lag analys is from giving defin i- 
tive results on its possible short nature (Norris et al. 2010b). 
The b urst is however considered a SGRB in iNorris et al. 
l|201ll ). 



2.1 Swift-BAT data analysis 

BAT data have been processed using standard Swift-BAT 
analysis tools within heasoft (v. 6.10). In particular, the 
BATGRBPRODUCT script has been used to generate event lists 
and quality maps necessary to construct 4 ms mask- weighted 
and background-subtracted light-curves in the 50-100 keV 
and 100-200 keV anergy bands. The ground-refined coor- 
dinates provided by the BAT-refined circulars have been 
adopted; standard filtering and screening criteria have been 
applied. 



2.2 Swift-XRT data analysis 

XRT data have been processed with the latest heasoft re- 
lease available at the time of writing (v. 6.10) and corre- 
sponding calibration files: standard filtering and screening 



■^ A 3(T threshold would only exclude GRB 051210, where the fluc- 
tuation has a significance of ^ 2.8a. 

^ The percentage of SGRBs with variable XRT light-curve 
8/60 ~ 13% is much less th an the ~ 30% of LGRBs showing 
flares jChincarini et al.ll2010l ). This result suggests that the per- 
centage of SGRB light-curves with variability superimposed is 
lower than in LGRBs. However, the lower statistics characteris- 
ing the SGRB curves prevents us from drawing firm conclusions. 
This topic will be addressed in a separate work. 
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Figure 1. Upper panel: 0.3-10 keV count-rate light-curve of 
GRB 100117A. Black solid line: continuous X-ray emission un- 
derlying the flare candidates computed as described in Section 
12.21 dashed lines: best-fitting flare candidate emission; red solid 
line: best estimate of the total emission. The vertical dot-dashed 
lines mark the flare candidate onset times. Inset: Complete Swift- 
XRT light-curve. The yellow flUed area marks the time window for 
the computation of the GCF lag (Sect. [Z3t . Middle panel: hard- 
ness ratio (HR) evolution with time; the HR is computed between 
1.5-10 keV (hard band) and 0.3-1.5 keV (soft band). Lower panel: 
Spectral photon index evolution with time as calculated by Evans 
et al., 2010. 



criteria have been appl ied. Pile-up corrections have been ap- 
plied when necessary l| Romano et al.l 120061 : IVaughan et al.l 
120061 ). Count-rate light-curves have been extracted in the 
total XRT 0.3-10 keV energy band as well as in the 0.3-1 
keV, 3-10 keV, 0.3-1.5 keV, 1.5-10 keV and 4-10 keV en- 
ergy bands. The 0.3-10 keV count-rate light-curves have 
been re-binned at a minimum signal-to-noise ratio SN=4 and 
then searched for statistically significant temporal variabil- 
ity superimposed over a smooth afterglow decay. A two-step 
procedure has been followed: first the smooth continuum 
contribution ha s been determined applying the method by 
iMargutti et al.l ([20111). A simple power-law or a smoothly 
joined broken power-law model is adopted (black solid line 
of Fig. [!}. As a second step, the properties of statistically 
significant fiuctuations with respect to the continuum have 
been determined adding a number of INorris et al.l (|2005l ) 
profiles to the best fitting continuum model. The best fitting 
INorris et al.l l|2005l 'l profiles constitute the sample of X-ray 
flare candidates of SGRBs analysed in this work. Figure [l] 
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shows GRB 100117A as an example: 3 distinct episodes of 
variability have been identified. The best fitting parameters 
of the entire sampl e are listed in Table [All The choice of the 
iNorris et alj (|2005l ) profile allows us to perform a one-to-one 
comparison with t he properties of X-ray flares and prompt 
pulse s in LGRBs (jChincarini et alJ I2OI0I : iBernardini et alJ 
l201lf ): Fig. [2] shows the evolution of the SGRB flare candi- 
dates width with time compared to LGRB flares. 

The evolution of the spectral properties of any source 
can be constrained through the analysis of its hardness ratio 
(HR), which is here defined as HR = clZl%'l-f''X^v) ■ ^ 
different binning with respect to the total 0.3-10 keV light- 
curve has been used for the 1.5-10 keV and 0.3-1.5 keV 
light-curves to improve the HR signal-to-noise ratio. The 
temporal evol ution of the s pectral photon index V has been 
calculated by Ev ans Il20ld . Results are portrayed in Fig. [U 
[BTl[B2llB3l andlB4l 

Count-rate light-curves have been converted into flux 
and luminosity (when possible) curves using the spectral 
information derived from a time-resolved spectral analysis 
where the spectral evolution of the source, if present, is prop- 
erly accounted for. This procedure allows us to convert the 
best fitting peak count rates of the X-ray flare candidates 
[A parameter of Table IA1|) into 0.3-10 keV peak luminosi- 
ties Lpk,?^ when the redshift of the source is known. Figure 
[3] shows the evolution of the SGRB flare candidate I/pk,F^ 
wi th time compared to the results obtained for LGRB flares 
bv IChincarini et al.l (|2010l ): a comparison of the two distri- 
butions can be found in Fig. [l] The isotropic equivalent en- 
ergy _Ejso.f^ h as be en determined integrating the best-fitting 
iNorris et al.l (|2005l ) luminosity-calibrated profiles from the 
onset time (fs) to ta + IQQw (where w is the flare candidate 
width). The uncertainty arising from the spectral calibration 
has been propagated following standard practice into the fi- 
nal I/pk,™ and Ef^jP u ncertainties listed in Ta bles [1] and 
W\\ We refer the reader to lMargutti et al.l (|2010al ') for details 
on the light-curves and spectra extraction. Figure [S] shows 
the SGRB flare candidates Ef^^ distribution compared 
to th e values dctcrnuned for LGRB flares, as computed by 
IChincarini ot al., (2010,). 



2.3 Spectral time-lag computation 

The spectral lag is the time difference between the arrival 
of high-energy and low-energy photons. For each GRB, the 
X-ray and 7-ray spectral lags (rj^g and t?, respectively) 
and associated uncertainties have been determined using a 
cross-correlation function (CCF) analysis. The CCF analysis 
requires the observations to have a fractional exposure equal 
to 1: this requirement excludes most of XRT observations 
taken in photon counting (PC) modqj. Among these, the 
late-time (f ~ 5 x lO" s) re-brightening of GRB 050724 (Fig 
IB1|) . We closely follow the prescriptions bv lStamatikos et al.l 



* Si«J/i-XRT automatically switches to the PC observing mode 
for count-rates below a few count s~^ to minimise the presence of 
pile-up. In PC mode, it is not uncommon to have short time inter- 
vals of no observation even during a single orbit. While the light- 
curve and spectra extraction procedures are basically insensitive 
to these short pauses, the CCF analysis would give un-reliable 
results. 
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Figure 2. Rest frame width vs. peak time relation for LGRB 
early-time flares from Chincarini et al., 2010 (open triangles) 
and SGRB flare candidates with and without extended emis- 
sion (red open and filled circles, respectively) . Inset: Complete 
view of the «'/(! -f z) vs. tpk/(l + z) relation established by 
LGRB flares obtained joining the data from Chincarini et al., 
2010 and Bernardini et al., 2011. The blue dashed line in both 
plots marks the best fitting relation calculated on LGRB flares: 

.J£- = 10(-i-0±0-5)('iEiLV^'^'^°'^\ where w and tpk are mea- 
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Figure 3. 0.3-10 keV peak luminosity evolution with time 
for LGRB flares (black open triangles, from Chincarini et al. 
2010) and SGRB flare candidates with and without extended 
emission (red open and flUed circles, respectively). Blue dashed 
line: best fltting power-law model for LGRB flares: -L'"^™ = 



]^q54.8±0.4 / 'pk 
' 1 + z 



and extrinsic scatter a = 0.73 it 0.08. 



Orange dot-dashed line: best fltting ip^ f^ decay re-normalised 
by a factor 100 to match the observed SGRB flare candidates 
rSGRB 
pk,F ■ 



(|2009l ) and lUkwatta et al.l (|2010l ') for the CCF computation: 
in particular, each CCF peak has been fltted using a third 
order polynomial; the number of points to be fitted around 
the CCF peak has been allowed to vary from case to case 
with the possibility to specify asymmetric intervals around 
the peak. In our analysis, a positive spectral lag is obtained 
if high energy photons lead low energy photons. 

The lag extraction is sensitive to a number of parame- 
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Figure 4. 0.3-10 keV (observer frame) isotropic equivalent peak 
luminosity ip^p^ of LGRB flares from Chincarini et al., 2010 
(hatched histogram) compared to SGRB flare candidates (fllled 
histogram); two vertical arrows mark the position of flare can- 
didates in SGRBs with extended emission. The vertical dashed 
lines mark the median values of the two distributions: ^p^p^ ~ 
lO'^Olergs-l; L^PS^ ~ lO'*'^-* ergs"!. Inset: 1-10000 keV rest 



ovrT o^l ■ T SGR.B 

frame isotropic equivalent peak luminosity distribution of the long 
(Nava et al., 2008) and short GRBs prompt emission (Ghirlanda 
et al., 2009; 2010), with median value: Lpk,p ~ lO^^-S ergs"!. 
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Figure 5. 0.3-10 keV (observer frame) isotropic equivalent energy 
EJ^GRB of LGRB flares from Chincarini et al., 2010 (hatched his- 
togram) compared to SGRB flare candidates (fllled histogram); 
two vertical arrows mark the position of flare candidates in 
SGRBs with extended emission. Vertical dashed lines: median 
EjLGRB ^ ]^q51 gj.g a^nj SjS^fB r^ lO^'^erg values. Inset: 1-10000 
keV rest frame isotropic equivalent energy distribution for the 
prompt emission of the widest samples of long (hatched his- 
togram) and short GRBs (filled histogram) with flrm spectral pa- 
rameter estimates at the time of writing (Amati et al., in prep.). 
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ters: energy band pass of each comparative light-curve, tem- 
poral bin resolution, signal-to-noise ratio and presence of 
background emission (i.e. in X-rays, the smooth X-ray decay 
underlying the time- variable signal). For the prompt 7-ray 
phase, the lag has been calculated using 4 ms light-curves 
(Section [23]) in the 50-100 keV and 100-200 keV energy 



bands. Time intervals covered by extended emission have 
been excluded. This allovifs us to perform a direct compari- 
son with the time-lag va lues obtained for LGRBs observed 
by BAT (Ukwatta et al.|[2010, ). Results are listed in Tabled 
column 5: the 8 SGRBs exhibit negligible r7 . 

In the X-rays the situation is complicated by the pres- 
ence of a smoothly declining afterglow emission underly- 
ing the episodes of possible activity (see e.g. black solid 
line of Fig. [T]). Choices of re-binning time-scales, energy 
bands and/or temporal intervals giving origin to correla- 
tion values (GCF peak) < 0.4 have been discarded. The 
choice of the energy bands to be compared is limited by 
the XRT 0.3 - lOkeV coverage. For each SGRB, the X-ray 
time lag ri^g has been computed for different energy bands, 
giving consistent results: the 0.3 — 1 keV and 3 — lOkeV en- 
ergy bands have been finally chosen to perfo rm a one -to 
-one co mparison to the results obtained by iMargutti et al.l 
(|2010bl ) for flares det ected in LGRBs . To this end, the LGRB 
flare time lags from IMargutti et al.l l|2010bl ) have been re- 
calcu lated using the GCF a, nalysis above (black dots in Fig. 
\E^: in lMargutti et al.ll2010bl a pulse peak lag was instead cal- 
culated (grey dots in Fig. ^. The pulse peak lag is defined 
as Ti^g'''' = ipcak - ipeak whcrc t^^^^ aud tf^^i, are the peak 
times of the best fitting profiles in the energy bands I and 
II, respectively. As such, t[^J^ is sensitive to the assumed 
pulse fitting model: while the dependency is limited in cases 
of bright events, the limited statistics of the SGRB X-ray 
light-curves would cause the pulse peak lag computation to 
be inaccurate. For this reason we refer to the CCF time lag 
for both short and long GRB data, in the 7-ray and X-ray 
regimes. The light-curve time binning can potentially affect 
the derived t^^^: for each SGRB the lag has been computed 
on light-curve pairs with 10 different time binnings spanning 
the range 0.2 — 20 s. The optimal time binning is defined as 
the lowest time scale giving origin to a GGF peak > 0.4 and 
is listed in Table [T] Larger binning time scales have been 
checked to produce consistent lag results. The window of 
time of investigation (ti and f/ of Table [l]) has been de- 
termined selecting the time interval containing positive, at 
least 1 — a signi ficant fluctuations ar ound the smooth X-ray 
continuum (see IMargutti et al.ll201ll for details). 

For each SGRB, ti and f/ have been varied of ~ 20% 
both towards larger and lower values: consistent time lag 
values have been found. The sensitivity of the lag measure- 
ment to the smooth X-ray light-curve decay underlying the 
candidate flares has been investigated calculating the lag 
''"lag" '-"^ light-curve pairs where the contribution of the 
smooth afterglow component has been properly subtracted 
and uncertainti e s prop agated following the prescriptions by 
IMargutti et al.l l|201ll ). For each SGRB this procedure has 
led to consistent Tj^g" — Tj^ values (Tj^g" systematically 
has larger uncertainties due to the lower signal-to-noise of 
the subtracted light-curves). For this reason we refer to rj^g 
hereafter. Finally, we have tested and verified the robustness 
of our choice of energy bands to be compared, re-binning 
times and window of time of investigation, by performing a 
number of simulations where artificial lags have been first 
introduced into the light curves and then successfully recov- 
ered. Results are reported in Table[T]. Figure |6]combines the 
Ti^g and Lpit luminosity information in the lag-luminosity 
plane and clearly shows that SGRB flare candidates fall off 
the prediction based on LGRB flares and prompt pulses. 
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Figure 6. Lag-Luminosity plot. Red circles: CCF lag for candidate flares of SGRBs: open symbols refer to SGRBs with detected EE; 
black (grey) triangles: CCF (pulse peak) lag for the sample of 9 flares of long GRBs of Margutti et al., 2010b, their Fig. 15. Black stars: 
prompt 7-ray data from the gold and silver sample of Ukwatta ct al., 2010. Red squares: 3(7 upper limits to the prompt lag of SGRBs 
for which it is possible to estimate the peak luminosity: open symbols refer to SGRBs with detected EE. The isotropic peak luminosity 
is computed in the 1 — 10^ keV and 0.3 — lOkeV bands for prompt data and X-ray flares, respectively; the lag corresponds to the time 
difference between light-curve structures in the 50-100 keV and 100-200 keV channels (prompt data) and 0.3-1 keV and 3-10 keV channels 
(X-ray flares). Blue dashed line: best fitting law for the LGRB prompt emission data. 

Table 1. X-ray and 7-ray properties for the sample of SGRBs analysed in this work. From left to right: GR,B name: a (*) indicates an 
early-type host galaxy morphology (Fong et al., 2010; 2010b), while GRBs with detected extended emission are in boldface; redshift, 
duration and average spectral photon index of the prompt 15-150 keV emission as determined from GCNs; extended emission (EE) 
duration and EE to IPC (Initial Pulse Complex) intensity ratio (i?int) from Norris et al. 2010a: upper limits on -Rint are listed when 
no EE has been detected (EE with zero duration); t; and tj define the temporal window for the X-ray lag calculation while the Atj^b 
column reports the time scale used to re-bin the X-ray light-curve pairs; CCF time lag computed between 0.3 — 1 keV and 3 — lOkeV; 
0.3 — 10 keV isotropic equivalent peak luminosity in the time interval ti — tf as determined from the Norris 2005 profile fit (see Table 
|AT}; short lived (SL) or long lived (LL) X-ray afterglow according to the classification by Sakamoto & Gehrels, 2009. 
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3 RESULTS 

The data analysis of the previous sections leads to the fol- 
lowing results: 



• SGRB flare candidates appear both in early-type and 
late- type host galaxy environments, irrespective of the short- 
lived (SL) or long-lived (LL) nature of the X-ray afterglow. 

• Both SGRBs with a bright extended emission (EE) and 
SGRBs which lack this component show cases of statistically 



significant fluctuations superimposed over smoothly decay- 
ing X-ray light-curves (Table [l|. 

• Flares in LGRBs are known to show a spectral hard- 
ening during the rise time and a softening during the 
decay time, reminisc ent of the prompt emission (e.g. 
iMargutti et al.l l2010bl and references therein) : as a result 
the hard ness ratio (HR) e volution mimics the flare profile 
(see e.g. iGoad et al.l 120071 . their Fig. 9) while the spectral 
photon index evolution anti-correlates with the flare flux. In 
spite of the lower statistics and limited AF/F of SGRB flare 
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Figure 7. Panel (a): Relative variability flux AF/F vs. relative variability time-scale At/t = wi/tpk for the sample of X-ray flare 
candidates in SGRBs (flUed circles), compared to early and late time LGRB X-ray flares (blue open diamonds and light-blue stars, 
respectively) from Chincarini et al., 2010 and Bernardini et al., 2011. The late time re-brightening detected in GRB 050724 is also shown 
for completeness with an orange open triangle (Bernardini et al., 2011). A small black dot marks data coming from SGRBs with detected 
extended emission. Solid, dashed and dot-dashed lines mark the kinematically allowed regions in different scenarios according to loka et 
al., 2005, their equations (7) and (A2). The At/t and AF/F distributions are portrayed in panels (b) and (c) adopting the same colour 
coding. The green tick marks in panel (c) show the flux contrast for the sub-sample of LGRB flares of Fig. |6] 



candidates (Fig. (Tjl , we find in the case with the best statis- 
tics a hint for correlation between the HR evolution and 
the temporal profile of GRB 101117A, with a photon index 
evolution that anti-correlates with the flux of the flare can- 
didates (Fig.[T] middle and lower panels). In the other cases 
(Fig. |BT] |B2] |B3l IB4|I the hmited statistics prevents us from 
drawing firm conclusions. 

• The SGRB fiare candidates width evolution is roughly 
linear in time and consistent with the w/{l-\-z) vs. ipk/(l+2:) 
relation established by LGRB fiares over ~ 4 decades in 
time (Fig. \^. The best-fitting law reads; w/(l -\- z) = 



10" 



'(tpk/(i + z)Y 



It is remarkable that data 



coming from LGRB flares both at early and very late time 
(beyond ipk/(l + z) ~ 10^ s, Fig. [51 inset) as well as tempo- 
ral fluctuations in completely different systems like SGRB 
afterglows are consistent at zero-or der with the same, ap- 
proximately linear, law. We refer to [Bernardini et al.l l201ll 
for a discussion of the possible bias affecting the w vs. tpk 
relation. 

• SGRB flare candidates are ~ 100 times dimmer than 



LGRB flares at the same rest frame time (Fig. [Sjl . Selecting 
the sub-sample of LGRB flares detected in the same rest 
frame time interval 60 s < tpk/(l + z) ^ 250 s of SGRB flare 
candidates, we obtain a median (I/p^p^) ~ 10'*^'*ergs~'^ to 
be compared to (Lpk,™) ~ lO^'^'^ergs"^ of the SGRB sam- 
ple showed in Fig. U) As a result, SGRB flare candidates 
fall off of a factor ~ 100 the peak luminosity vs. time re- 
lation established by LGRB flares which reads: -Lpk f^ ~ 

1054'«±"-4(ipk/(l + ^))"'''*°-'- 

• Short and long GRBs show a comparable 1 — 10* keV 
(rest frame) isotropic peak luminosity during their prompt 
emission, with a median (Lpk,p) ~ lO®'^'^ ergs~^ (Fig. |H in- 
set). On the contrary, the peak luminosity of flares of both 
categories evaluated at the same rest frame time 60s < 
tpk/(l + z) < 200 s differ of a factor ~ 100 as noted above. 
While for LGRBs the typical prompt (1-10'' keV) to flare 
(0.3-10 keV) peak luminosity ratio (Lp^P^/ip^p'^) ~ 300, 
for SGRBs the same quantity reads: (I/pk^p Apk,™) ~ 
3 X 10*. Flare candidates in SGRBs are therefore less lu- 
minous than expected using the prompt-to-flare luminos- 
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ity scaling observed in LGRBs at the same flare rest frame 
tims0. 

• SGRB flare candidates are ~ 100 times less energetic 
than LGRB flares (Fig. [SJ, with a median 0.3-10 keV en- 



Table 2. Summary of the properties of SGRB flare candidates 
compared to LGRB X-ray flares observed at the same t-py^/(l + z). 



ergy (i?S?|«) ~ 10«-« erg ((£L"^"> ~ lO^"^'^ erg). Since the 
width of LGRB flares and SGRB flare candidates are com- 
parable when evaluated at the same tpk/(l + z), this result 
is a natural consequence of the (Lp°™)/(ip°F^> ~ 100 
On average, in the time interval 60s < 



reported above. 

tpk/(l + z) < 250 s a flare in a LGRB emits 



0.6% of the 



1 - 10" keV prompt BisS^p , the observed SGRB flare can- 
didates isotropic energy is ~ 0.2% the 1 — 10* keV prompt 

TTiSGRB |6l 

• Flares and prompt pulses in LGRBs define a spectral 
lag-peak luminosity relation (Fig. [6l): this finding is high ly 
suggestive of a common origin (jMargutti et al.l l2010bl rl. 
On the contrary, SGRB prompt pulses are known to ex- 
hibit much shorter lags than expected if they were to fol- 
low the LGRB pro mpt pulses lag-luminosity relation (e.g. 
iGehrels et al.ll2006h . Figure[S]extends this behaviour to their 
flare candidates: like SGRB prompt pulses, flare candidates 
in SGRBs fall off the lag-luminosity relation defined by 
LGRBs. 

• Flares candidates of SGRBs in the AF/F vs. rela- 
tive variability time-scale (Ai/i = w/tpk) plane are com- 
patible with variability arising from density fluctuations of 
many regions viewed off-axis: on the contrary, neither the 
refreshed-shock nor the patchy-shell scenario is able to ac- 
count for the ob served properties of the entire sample (see 
lloka et al.J|2005l for details on the deflnition of the various 



3I. 



scenanos]J. In particular: a K-S test comparing the At/t 
distributions of LGRB flares and SGRB flare candidates re- 
veals that they belong to the same parent population at 
~ 10% level of probability; the probability reaches the 88% 
level if LGRB flares are selected in the SGRB peak time 
range (60s < tpk/(l + z) < 250s). This is consistent with 
the common ui/(l + z) vs. tpk/(l + z) relation followed by 
short and long GRBs discussed above. On the contrary, no 



® The prompt emission peak luminosity is likely to be biassed 
towards the bright end of the real Lpk,p distribution, since a min- 
imum signal-to-noise is needed to constrain the spectral parame- 
ters and calculate Lpk,p in the 1-10* keV range: this requirement 
is more severe in the case of SGRBs whose observed emission is 
usually less bright than LGRBs. However, Fig. [8] shows that the 
main conclusion of this paragraph remains true even after relax- 
ing the requirement above. 

® Again, this ratio is likely to be a lower limit to real value due 
to the bias affecting the the sample of SGRBs with prompt E\^a 
measure discussed in the previous paragraph. 
^ Note that flare and prompt lags are calculated in different en- 
ergy bands and are not directly comparable. The fact that flares 
define a lag-luminosity relation with slope very similar to the 
prompt data is however sugge stive of a strict connection between 
flares and prompt pulses. See lMareutti et al.ll2010b for details. 
* The smoking gun against a refreshed shock scenario would be 
the detection of a spectral change contemporaneous to the flare 
candidates: while the statistics of the XRT light-curves of SGRBs 
is limited, in the case with best statistics the HR evolution is 
correlated to the candidate flare profiles and the continuum after 
the flaring emission is softer than the emission detected during 
the period of temporal variability (Fig. [TJ . These findings favour 
an alternative explanation. 
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X-ray flare candidate in a SGRB shows a relative variability 
flux AF/F > 2 in strong contrast with the LGRB AF/F 
distribution at comparable ipk/(l-|-z) (Fig.[7ll. Such promi- 
nent flares would be easier to detect, so that it is unlikely 
that an observational bias could explain the present lack 
of detection. A K-S test comparing the two distributions 
shows that the probability that LGRB and SGRB flare can- 
didates share the same AF/F parent population is as low as 
3.3 X 10"'*. This result partially inherits the uncertainty af- 
fecting the completeness of both distributions for very small 
AF/F values. Another source of uncertainty arises from the 
difficulty in evaluating the continuum underlying the pos- 
sible flare emission in SGRBs when data are particularly 
sparse. In spite of these limitations, after more than 6 years 
of Swift observations (and ~ 60 SGRB afterglows detected) 
there is still no SGRB showing a prominent {AF/F > 10) 
fast-varaibility At/t ^ 1 feature during its X-ray afterglow. 
The SGRB flare candidate AF/F is instead more similar to 
the relative variability flux of flares in LGRBs detected at 
late times (tpk > 1 ks, light-blue stars in Fig. [6l main panel; 
iBernardini et al.ll2011l ): the two AF/F distributions share 
the same parent distribution at ~ 21% level of probability 
(K-S test). 

The results above demonstrate the complexity charac- 
terising the SGRB flare candidates phenomenology: Table 
[2] reports a summary of their properties when compared 
to LGRB fiares. The main result is that the population of 
SGRB X-ray fiare candidates only partially share the set 
of observational properties of LGRB X-ray fiares detected 
at the same rest frame time: are there real X-ray flares in 
SGRBs? A detailed discussion is provided below. 



4 DISCUSSION 

Observations show that like LGRBs, at least some SGRB X- 
ray afterglows deviate from a smooth power-law decay and 
show variability. In the following we discuss the properties 
of SGRB fiare candidates providing a one-to-one comparison 
with LGRB X-ray fiares. The aim is to better understand the 
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origin of the short burst afterglow variabihty and uncover 
potential links with the prompt phase. 

With Af/i > 1, the late-time (ipk ~ 5 x 10* s) 
re-brightening of GRB 050724 strongly differs from the 
properties of the entire sample of flare candidates de- 
tected in SGRBs thus questioning its classification 



flare-like ep i sode ( Panaitescu et al. 200d : see however 



Grupe et all l2006l: ICampana et alJ l2006l : iMalesani et all 



upe -.-y 

2007l 'l. iMalesani et al.1 (|2007l ') report the detection of an 
optical and radio re-brightning associated to the X-ray 
bump which is unusual if compared to the standard prop- 
erties of X-ray flares, while being more common to late- 
time re-brigh tneings observed in L GRBs as well (see e.g. 
GRB 081028. iMargutti et al.ll2010al '). In addition, n o hard 



to-so ft evolution can be detected in the X-ray data (lEvans I 
bold ), which is instead typical of flares ("Goad et al.1 120071 ) 
and prompt pulses in LGRBs (Hakkila ct al. 2011). In the 
following we focus our attention on SG RBs fast variabil- 
ity (A t/t < 1) referring the reader to iBernardini et al.l 
(|201ll ) for a complete discussion of the late-time behaviour 
of GRB 050724. 



4.1 SGRB local and global environment 

The standard model fsee lNakar et al.ll2007l for a recent re- 
view) explains the X-ray afterglow of long and short GRBs 
as synchrotron radiation arising from the deceleration of 
a relativistic blast wave into the external medium. If the 
shock front is homogeneous and expands into a smooth am- 
bient density, a smooth afterglow light-curve is expected. 
In this context, variability in the X-ray afterglow could be 
caused by re-freshed shocks (i.e. shocks caused by slow shells 
catching up with the leading, decelerated shell at late times, 
iKumar fc Piranll2000al . iGranot al.ll2003l ): Fig.[7]shows that 
half of the SGRB flare candidates sample do not comply with 
this scenariqj. Furthermore, the spectral variability shown 
in Fig. [ll lower panel, makes it difficult to interpret the flare 
candidates in the re-freshed shocks scenario. 

A first alternative is to re lax the assumption on the 
homogeneity of the shock front (|Kumar fc Piranll2000bl ) : an 
intrinsic angular structure of the emitting surface is able to 
produce variability with a characteristic time-scale At ^ t 
if th e angular fluctuation is persistent (patchy shell model, 
iNaka r ct al. 2004l^"l: no SGRB flare candidate is consistent 
with this expectation (Fig. [T] main panel). 

A second alternative invokes the presence of ambient 
density fluctuations either caused by turbulence in the ISM 
or by variable winds from the progenitor. From Table [l] 
it is however apparent that temporal X-ray variability has 
been detected for SGRBs residing both in early-type and 
late-type host galaxies which likely have very different ISM 
properties. In particular GRB 050724 and GRB 1001 17A are 
the unique two SGRBs with secure early-type host galaxy 
association (JBarthelmv et al.ll2005bl : iFong et al.ll20ld ). This 



3 See however [Cranot al.l | |2003^ : Iviasis et al.l 1I2OIII) : sharp op- 
tical and radio flares could be produced by collision of ultra- 
relativistic shells. 

^^ Details on variability arising from a time-varyi ng anisotropic 
emitting surface can be found in llol<a et al.J 1120051 ). 



suggests that the ISM turbulence is unlikely to provide a fea- 
sible physical mechanism for the detected variability. Note 
however that the limited size of our sample prevents us 
from quantitatively discussing the correlation between the 
appearance of flare candidates and host environment. 

Different progenitor models of SGRBs lead to distinct 
predictions on their local environment as well. In particu- 
lar, according to the sta ndard compact binary merger sce- 
nario (NS-NS or NS-BH. lEichler et al.|[l989l : lNaravan et al.1 
[1993), no variable wind is expected from the progenitor. 
An origin of SGRB ffare candidates from density ffuctua- 
tions of the circumburst environment is therefore consid- 
ered unlikely. Alternative scenarios leading to SGRB envi- 
ronments with circumburst den sities comparable to LGR Bs 
have been however explored bv I Nvse wander et al.l (|2009l ) to 
explain the similar Fr/Fx ratio (where Fr and Fx stand 
for the afterglow flux density at 11 hrs post-trigger in the 
R and X-ray bands, respectively). A systematic difference 
between sub-galactic environments able or unable to pro- 
duce variability in the X-ray afterglow could in principle 
be revealed by different offsets with respect to the host 
galaxy centers: however, the present sample of SGRBs with 
flare candidates includes both SGRB wit h large offsets (e.g . 
6 = 14.80 ± 0.34 kpc for SGRB 071227, iFong et al.l [2010131 : 
iD'Avanzo et al.l |2009| ) and events with very small offset s 
(e.g. 5 = 0.47 ±0.31 kpc, SGRB 100117A. lFong et allboid ). 
While the observed offset distribution is incomplete, these 
data suggest that the properties of the local environment of 
SGRBs are not the key parameters determining the presence 
of flare candidates in their X-ray light-curves. This conclu- 
sion is strengthened by the results of Sect. 14.51 

Flares candidates and prompt pulses could alternatively 
share a common origin. In this case, flare candidates would 
bring no information on the external medium density. We re- 
fer to this interpretation as the "internal origin" possibility 
(see Sect. Ft)) . 

4.2 Flare candidates and extended emission 

3 bursts in our sample (boldface in Table[T]) present evidence 
for extended emission (EE): a long-lasting {At ~ lO^s) soft 
X-ray tail that follows the short hard spike in the prompt 
phase. iNorris et al.l (|2010al ) analysed a sample of 51 Swift- 
BAT SGRBs looking for the presence of EE in their 7-ray 
data and concluded that ~ 3/4 of the BAT SGRBs are not 
accompanied by an EE component. In particular, in this 
work it is shown that the upper limit on the EE to IPC 
(initial pulse complex) intensity ratio of the SGRBs of our 
sample (i.e. with flare candidates) not showing EE is a factor 
> 10 below the sta ndard values found for GRBs with EE 
(jNorris et al.ll2010al . their Fig. 1). This implies that in those 
cases the EE is either very weak or absent. The inverse is 
also true: some SGRBs with bright EE are accompanied by a 
smooth X-ray light -curve at t > 80 s (see e. g. GRB 080503, 
iNorris et ai]i2010al . their Table 1). 

While it is still unclear if the origin of the EE compo- 
nent is related to the prompt emission, the afterglo w or a 
third, unknown physical process (Perlcy 2009; Norr is et al.l 
l2010al : INorris et al.ll201ll ). our analysis shows that the pres- 
ence of a bright EE does not imply the presence of flare 
candidates and viceversa. We stress that the higher aver- 
age brightness characterising the XRT light-curves of GRBs 
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with EE with respect to those without EE, would naturally 
bias the result in the opposite direction (i.e. it would favour 
temporal fluctuations to be detected in XRT light-curves 
with the best statistics, leading to a biassed flare-EE con- 
nection) thus strengthening our conclusion. The limited size 
of our sample (which is however the widest possible at the 
time of writing) does not allow us to quantitatively discuss 
the flare candidates vs. EE correlation (or lack thereof). 

Our analysis cannot however exclude that SGRB flare 
candidates constitute temporal fluctuations superimposed 
on (and physically linked to) the X-ray tail of the EE, 
as suggested by the epoch of flare candidates detection 
(tpk/(l + z) < 200 s). Furthermore, the limited range 
spanned by the flare-to-continuum flux ratio AF/F (with 
AF/F ~ 1, Fig. [7| is suggestive of a physical link between 
flare candidates and the underlying continuum (i.e. the EE). 
Again, flare candidates would be associated to both bright 
and faint (undetected in the 7-rays) EEs. 



4.3 Time scales 

The relative time-scale distribution At/t of SGRB flare 
candidates (Fig. [T] upper panel) is compatible with being 
drawn from the same parent population of flares detected 
in LGRBs at similar tpk/(f + z) at ~ 88% level of proba- 
bility: with a median value of {At/t) = 0.25 and extending 
from O.f up to 0.5 it furthermore satisfies the limit At/t ^ f 
which is the classical argum ent against fluctuat ions arising 
from the external shock fe.g. lZhang et al. llJOO^ '). The At/t 
distribution therefore does not support an exter nal shock 
origin for flare candidates in SGRBs (see however iDermer I 
[2003), in agreement with the findings of Sec. 14.11 

The evolution of the flare duration with time «;(ipk) for 
tpk > T90 (Fig. [2]) observa tionally distinguishes fl ares from 
prompt pulses in LGRBs (jMargutti et al.ll2010a ): as time 
proceeds, LGRB flares becomes wider and wider following 
the best fltting law: w/{l -I- «) ~ (ipk/(l + z))^''^^^'^ ■ This 
quasi-linear regime establishes for t^y/{l-\- z) > r9o/(l + 2:), 
and extends 4 decades in time up to fpk/(l + 2) > 10^ s. No- 
tably, while the -Ei;J^p^ released during the LGRB prompt 
emission spans more than 3 decades in energy (Fig.[5|, likely 
reflecting different properties of the LGRB central engines, 
the duration of subsequent episodes of activity seems to fol- 
low a universal evolution with limited dependence on the 
initial energy budget. Flares with different amplitudes A, 
flux ratios AF/F and fluences show similar w/(l -\- z) whe n 
appearing at the same ipk/(l + z) (|Chincarim et al.ll2010h . 
suggesting that the main parameter driving the flare width 
evolution in a LGRB is the elapsed time from the explosion 

ipk/(l + ^0 

Figure O shows that flare candidates in SGRBs are con- 
sistent with the same quasi-linear temporal scaling: from the 
width measurement it is not possible to distinguish a flare 
in a LGRB from a flare candidate in a SGRB. The tem- 
poral properties of the prompt emission of long and short 
GRBs are clearly different in terms of duration and vari- 
ability (e.g. iNakar et al.ll2002l ): however, ~ 30 s later, the 



width of fluctuations superimposed over their X-ray light- 
curves seems to have lost any information on the nature of 
central engines able to produce 7-ray photons with such dif- 
ferent temporal properties. In both cases the w evolution is 
driven by the ipk, irrespective of their different initial con- 
ditions (and initial variability time scales): while for flares 
w|°"'^/(l + 2) ~ w'^'^^'^/{l + z) at similar ipk/(l + z), 
for prompt pulses wf^'^^ < w^'^^^ (with wf^^'^ ~ 0.05 
s to be compared to M;p ^^ '7:' ^ ^' ob server frame values, 
iNakar et al.ll2002l : iNakar et al.ll2002tf ^1 ). 

This flnding strongly suggests that the origin of the 
quasi-linear w/{l + z) vs. fpk/(l-f2:) scaling must be within 
what is in common for the long and short GRB models, 
irrespective of the progenitors, environment, lifetime, vari- 
ability time-scales and prompt energy release of their cen- 
tral engines. Both models require the prompt emission to 
originate from a relativistic outflow (see e.g. |Piran 2004J): 
if the longer and longer duration of flares is of dynamical 
origin and dominated by the expansion of the emitting re- 
gions, no memory of the properties of the central engine 
which launched the relativistic outflow would be preserved 
so that long and short GRBs would display flares with sim- 
ilar u;(tpk). Alternatively, in the context of accretion mod- 
els the Ui(tpk) relation originates from the viscous evolution 
of the hyper-accreting disk around the newly born black 
hole, co mmon i ngredient of the likely progenitors of the two 
classes (jPerna et al. 2006 ). Our results would imply a strict 
analogy between the mechanisms that regulate the late-time 
evolution of accreting disks originating from coUapsars and 
mergers of compact objects, irrespective of their different 
masses and life times (according to the standard scenario 
accreting disks related to SGRBs are likely to be less mas- 
sive then LGRB disks and sh ort-lived, as sugge sted by the 
obs erved prompt duratio n. See lNakar et al.ll2007l for details). 
iLazzati et al.l ([20111) recently suggested that instabili- 
ties arising from the propagation of the jet through the dis- 
rupting star could explain the presence of flares in LGRBs 
with AF/F < 10. Even assuming a continuous and feature- 
less release of energy by the central engine, the pressure of 
the surrounding stellar material would provide the physical 
origin for the intermittent flare behaviour, naturally explain- 
ing the universal (i.e. with limited sensitivity on the star 
roperties and ener gy budget) quasi-linear «;(fpk) relation 
Lazzati et al.ll201ll . their Eq. (6)). However this model can- 
not account neither for presence of flares in SGRBs nor for 
their similarity to LGRB flares in the ■w/{l-\-z) vs. fpk/(l+2) 
plane: according to the majority of SGRB progenitor scenar- 
ios the engine is "exposed" and there is no stellar material 
the jet could interact with. As a consequence, if the u;(ipk) 
relation in LGRBs originates from the interaction with the 
progenitor stellar material, it is difficult to explain why flare 
candidates in SGRBs, while originating from a completely 
different physical scenario, are however consistent with the 
same relation, as observed. Our results therefore imply that 
either the LGRB flare «;(ipk) relation does not originate 
from the jet-stellar material interaction or that the progen- 



^^ As a second level of approximation, one should cons ider that 
more p rominent flares (A ^ 1) tend to be wider (.Margutti et al.l 
l2010bl) . 



^^ Note however that the first 1-2 s of bright LGRBs display 
variability time sca les comparable to the SGRB prompt emission. 



iNakar et afl ||2002| '|. 
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itors of long and short GR Bs are much more similar than 
previously thought (see e.g. iLazzati et al.lboid 'l. 

Finally, our findings suggest that while the variability 
time scales measured during the prompt emission could still 
dire ctly refiect t he intrinsic variability of the central engine 
(see[Piran|[200J and references therein), absolute measures 
(i.e. not re-normalised) of flare widths likely do not (if time- 
dilated by physical mechanisms which are only indirectly 
related to the central engine). On the contrary, the ratio of 
interesting time-scales of the same flare profile, being sub- 
ject to the same temporal stretching, could still retain an 
imprint of the original mechanism at work: this would ex- 
plain why in LGRBs the flare rise ti me tr and decay time td , 
like the w, linearly ev olve with fpk (jChincarini et al.ll2010l : 
iBernardini et al.ll201ll ) preserving their ratio U ~ 0.5td over 
4 decades in ipk and leading to flares with asymmetry val- 
ues very similar to the prompt pulses (while being a factor 
> 100 widerQ 



4.4 Flux contrast 

While the «;(ipk) and At/t measurements do not allow us 
to distinguish a flare in a LGRB from a flare candidate 
in a SGRIO, the flux contrast distributions AF/F of the 
two populations are significantly different (Fig. [7|), with 
SGRB fiare candidates having systematically lower AF/F 
values. In contrast to LGRBs, none of the ~ 60 Swift X- 
ray afterglows associated to SGRBs shows cases of powerful 
{AF/F > 5), highly variable At/t < 1 flares. A Chandra 
observation of 9 X-ray photo ns from SGRB 5070 JHI ~ 15 
days after the explosion led I Fox et al.l ((20051) to conclude 
the presence of high-amplitude {AF/F « 10), fast variabil- 
ity {At/t ~ 0.01) in its X-ray afterglow. Swift observations 
suggest that this kind of variability, if present, is extremely 
rare. SGRBs lack the presence of prominent fast-rise expo- 
nential decay features superimposed over their X-ray after- 
glow foi_tgk/(]_+2:)^60 s. 

iLazzati et al.l (|2011l ) predicted that if flares in SGRBs 
originate from the intrinsic variability of their inner engine, 
their AF/F distribution should be analogous to the brighter 
population {AF/F ~ 100) of LGRB fiares. With a maxi- 
mum AF/F ~ 2 the detected SGRB fiare candidates are 
not consistent with these expectations and populate instead 
the low end of the LGRB distribution unless the SGRB con- 
tinuum flux F has been overestimated by a factor ~ 100 
which we consider unlikely. An interesting possibility is that 
the X-ray light-curves of SGRBs are dominated by an emis- 
sion component which is not present in the LGRB after- 
glows (such as the EE): this would lead to systematically 
lower AF/F for SGRBs when compared to LGRB fiares. 



^^ Note that given the limited statisties of the SGRB X-ray af- 
terglows, nothing can be said about the asymmetry of SGRB flare 
candidates. 

^* The limited number of flare candidates plays a role in this 
statement. We cannot exclude that a significant improvement of 
the SGRB and LGRB statistics could lead to the detection of 
secondary effects. 

1^ GRB 050709 is a short burst detected by HETE for which 
Swift did not do the follow up. For this reason it is not included 
in the present sample. The episode of prominent variability is 
possibly related to a statistical fluctuation. 



A LGRB f lores 1|„<1 ks 

X LGRB f lores t|„>1 ks 
• SGRB flare condidotes 
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Figure 9. Relative variability flux AF/F evolution as a func- 
tion of tpk/Tgo: black triangles and green stars indicate LGRB 
flares with ipk < 1 ks (Chincarini et al., 2010) and tpk > 1 ks 
(Bernardini et al., 2011), respectively. Red open and filled circles: 
flare candidates in SGRBs with and without extended emission, 
respectively. Horizontal dashed, dot and dot-dashed lines mark 
the AF/F median value for the three samples. 



The clustering of the AF/F ratio around 1 would suggest 
a link between fiare candidates and underlying continuum. 
However, the strong correlation found bv I Nvse wander et al.l 
( 2009J) between the X-ray afterglow fiux and prompt fiuence 
of both long and short GRBs reveals a high degree of similar- 
ity between their X-ray afterglows at least at late times (i.e. 
11 hrs, observer frame). A detailed comparison of short and 
long GRB X-ray afterglows at early times is in preparation. 
The SGRB AF/F observed at tpk ~ 100 s is instead 
typic al of LGRB fiares detec ted at much later times: tpk > 
Iks (Be rnardini et al.l 120111 ). From Fig. [5] it is apparent 
that SGRB fiare candidates are consistent with the AF/F 
distribution of LGRB fiares detected at the same tpk/Tgo 
epoch (where Tgo is the duration of the prompt 15-150 keV 
emission): differences instead arise if we compare the prop- 
erties of the two classes at the same ipk/(l + z). While 
the Tgo is possibly a crude estimation of the intrinsic time 
scale of evolution of the central engine, this result opens 
the possibility that prominent fiuctuations are not currently 
detected in SGRB afterglows because of the intrinsically 
late-time Swift re-pointings: trop ~ 100 s corresponds to 
trep/Tgo ~ 100 - 1000 (SGRBs) and to trep/Tgo - 1 - 10 
(LGRBs). From another perspective this finding directly 
connects the fiux properties of SGRB flare candidates to 
the evolution time scale of their central engines. The prompt 
Tgo qualifies as a good proxy for the intrinsic time scale that 
drives the subsequent flux evolution of flares and flare candi- 
dates in long and short GRBs, respectively. This conclusion 
is strengthened by the results of Sec. 14.51 



4.5 Energetics 

4-5.1 Flare Lpk evolution with time 

While SGRB prompt pulses compete with those of LGRBs 
in terms of peak luminosity, the same is not true for their 
late-time variability: SGRB flare candidates are a factor 
~ 10^ dimmer that expected (Fig. [Sjl. Sec. 14.41 and Fig. 
[5] suggest that long and short GRBs might be astrophysi- 
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Figure 8. Re-normalized peak luminosity vs. re-normalized peak time for the sample of LGRB flares from Chincarini et al. (2010) and 
Bernardini et al. (2011) (triangles) with tpk/Tgo < 10^ and SGRB flare candidates (circles). Dark and light colours distinguish events 
with and without redshift measurement to allow a direct comparison with Fig. |3] Open symbols refer to SGRBs with detected extended 
emission. Black dot dashed line: best fitting power-law model for the LGRB sample: Lpk,F/(ip) = lO""'^'^ (tpk/Tgo)""'® and 

extrinsic scatter a = 0.62 ± 0.01. The shaded area marks the ±1 a region around the best fit. From top to bottom, the blue dotted lines 
mark the i?iso,F/Siso,P = 1,0.1,0.01,0.001 regions of the plane, where i?iso,F and i?iso,P have been calculated in the 0.3-10 keV and 
15-150 keV energy bands, respectively. The flare width vs. peak time relation from Chincarini et al. (2010) has been used. 



cal systems evolving on completely different time-scales and 
with different energy budgets but otherwise based on a sim- 
ilar physical mechanism of emission. In that case, we would 
expect the evolution of the flare luminosity to exhibit much 
better agreement between long and short GRBs once the in- 
trinsic time-scale and central engine energy scaling are prop- 
erly accounted for. Figure [S] shows that this is indeed the 
case: long and short GRBs are consistent with a common (al- 
beit highly scattered) Lpk,F/(I/p) ~ (tpk/r9o)''''^='=°'^ scal- 
ing (where (Lp) and Tgo are the average prompt luminosity 
and duration of each GRB in the 15-150 keV energy range, 
respectively). Figure [8] therefore establishes a direct connec- 
tion between the properties of SGRB flare candidates and 
LGRB flares, providing further support to a common, likely 
internal origin (see Section [4. 7|) . 

In the context of a one to one comparison of short vs. 
long GRB X-ray variability, the re-normalised peak lumi- 
nosity and re-normalised peak time are the most relevant 
quantities: the observed Tgo and (Lp), with their sensitivity 
to the instrument threshold and energy band passes, likely 
provide crude (but nevertheless the best) approximations 
to the exact values to be used (while being partially re- 
sponsible for the large scatter of the relation). In particu- 
lar the loss of total fluence in the prompt 7-ray due to the 
limited Swift-BAT band-pass likely affects short more than 
long GRBs. As a result the SGRB I/pk,F/(T'p) values may be 
overestimated when compared to LGRB valuea^^l. However, 



This effects would mainly depend on the different prompt 



Fig. [5] shows that even a factor of several of overestimation 
(Nvs ewander et al.ll2009l ) would not strongly affect our con- 
clusions. 



4-5.2 Implications of the I/pk,F/(ip) vs. fpk/Tgo relation 

SGRB flare candidates are consistent with the liighly scat- 
tered Lpk,p/(//p) vs. ipk/Tgo relation established by LGRB 
flares: since the origin of LGRB flares is l ikely connected 
to the ir prompt pulses (|Krimm et al.l 120071 : iMargutti et al.l 
l2010bl ) , we speculate that a similar physical mechanism (ex- 
cept for its energy budget and life-time) powers long and 
short GRBs: observationally speaking, the main distinction 
between flares and prompt emission in LGRBs is the evo- 
lution of the former with time for ipk > Tgo- It is therefore 
natural to expect a similar scaling of LGRB flares and SGRB 
flare candidates in terms of tpk/Tgo if they share a common 
origin. 

A comparison of the prompt 7-ray emission of SGRBs 
to the initial 2 s of LGRBs reveals a high degree of sim- 
ilarity in the pulse duration distributions and correlations 
in the temporal structure of the two classes (JNakar et al.l 



emission spectral peak en e rgy E peak of long and short GRBs: 
however, iGhirlanda et al.l 1120091 ) showed that SGRBs have a 
harder low-energy spectral component but only slightly higher 
^poak when compared to LGRBs. 
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|2002| : iMcBreen et al.1 200lJM;_an analogous study was per- 
formed by iGhirlanda et al.1 l|2009f ) in the spectral domain: 
based on the spectral analysis of the prompt emission of 79 
short and 79 long GRBs detected by BATSE the authors 
showed that no difference is found comparing the spectral 
properties of SGRBs to the first 1-2 s emission of LGRBs. 
Temporal and spectral studies therefore point to a common 
mechanism operating in the first few seconds of any event. 
The present work extends this similarity to their late-time 
emission. 

The quasi-linear K;(tpk) evolution shared by long and 
short GRBs allows us to draw reference lines of equal flare- 
to-prompt energy ratios as a function of ipk/Tgo in the 
Lpk^F / {Lp) plane (blue dotted lines of Fig. [8]). Flares in the 
0.3-10 keV band pass are found to emit between (0.1 — 100)% 
of the prompt 15-150 keV _Biso,p, with the majority of them 
lying between the 1% and 10% levels. Little evolution of 
the flare-to-prompt energy ratio £iso,F/£'iso,p in terms of 
tpk/Tgo can be inferred from the plot (for fpk/Tgo < 300). 
In particular, SGRB flare candidates and LGRB flares show 
comparable flare-to-prompt energy ratios, as reported in Ta- 
ble[2l This finding provides further support to a physical link 
between LGRB flares and SGRB flare candidates. 

The GRB central engines seem to release comparable 
fractions of prompt emission energy at late times, irrespec- 
tive of the long or short GRB nature. From Fig. [S] 
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with a — —0.1 ± 0.1. The best-fitting rest-frame to(tpk) re- 
lation reads: 
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with /3 = 0.2 ± 0.2 (Fig. [2]). Equations (O and ^ express 
common properties of LGRB flares and SGRB flare candi- 
dates. The normalisation parameters A'^i and A'^2 possibly 
hide the dependence of I/pk,F/(ip) and w on other parame- 
ters. This hidden dependence might be partially responsible 
for the large scatter of relation ((l|; a and /? parametrise 
the non-linear dependence of I/pk,F/(ip) and w/{l + z) on 
tpk/Tgo and ipk/(l + z), respectively. Combining Eq. ([T]) and 
Eq. © it is possible to show that: 



Ei, 



Eis 



= M 



Ipk 



$-a 



1 + 2 



(3) 



wher e A/" ^ 0.9^jyiJV2 and Biso.F ~ 0.9LpkU)/(l + z) (valid 
for a iNorris et al.ll2005l flare profile where w is calculated be- 
tween 1/e intensity po ints and with tr = 0.5td as found by 
IChincarini et al.ll2010l for LGRB fiares). Equation ((3]) shows 



that Eiso,F / Eiso,p oc [fpk/(l-l-; 



1/3-a. 



a weak dependence of 



the flare-to-prompt energy ratio on ipk/(l + z) is expected 
if P — a r^ Q. From the best-fitting relations, we find that 
both parameters are consistent with at la: in particu- 
lar /3 — a = 0.3 ± 0.2 (consistent with at 1.5 (t level). At 
similar ipk/(l + z) a. residual dependence on the rest-frame 
prompt duration Tgo/(l + z) is expected to arise from the 

^^ The temporal evolution of pulses as a function of frequency 
(i.e. the spectral lag) shows however dissimilarities, as discussed 
in Section 14.61 



third term of Eq. ([3]): this dependence, if present, would 
be able to distinguish the population of LGRB flares from 
SGRB flare candidates detected at the same rest-frame peak 
time in terms of £'iso,F/£'iso,p. However, a = — O.liO.l (and 
the relation is highly dispersed). We therefore conclude that 
the quasi-linear tu(tpk) and I/pk,F/{ip) vs. tpk/Tgo relations 
translate into iSiso.F/S'iso.p ratios which, at flrst order ap- 
proximation, show limited dependence on both the proper- 
ties of the central engine (i.e. duration of the prompt emis- 
sion) and elapsed time from the explosion. 

From Fig.H Lpk.p'^/ip^F^ ~ 10"^ at the same tpk/(l+ 
2). A factor ~ 100 is roughly the ratio of the isotropic energy 
emitted by long and short GRBs during their prompt 7-ray 
emission (Fig. [5l inset J^. LGRB flares and SGRB flare can- 
didates are therefore expected to show a similar behaviour 
in the I/pk.p/fi'iso.p vs. tpk/(l + z) plane (in strict analogy 
with the X-ray afterglow scaling found by Nyscwandcr ct ajj 
|2009| . their Fig. 6). Equation 11} can be easily re-arranged 
into: 
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-(a + 1) 
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Again, the limited departure of I/pk,F/(ip) from a linear 
relation in ipk/Tgo (a = -0.1 + 0.1) causes LGRB flares and 
SGRB flare candidates to share the same scaling (at least 
at flrst order approximation) . 

All the above indications point to a SGRB flare candi- 
dates internal origin (Section l4.7|l and establish a connection 
between long and short GRB X-ray variability. 



4.6 The lag-luminosity relation 

Negligible spectral lag above ~ 25 keV is the fundamental 
attribute that makes the prompt 7-ray emission of short 
bursts different from LGRBs, in addition to their nar- 
rower pulses, shorter du ration and slightly harder emission 
(jNorris fc Bonnelll 120061 and references therein). The spec- 
tral lag is the time difference between the arrival of high- 
energy and low-energy photons: in our analysis, a positive 
value indicates that high energy photons lead the low en- 
ergy emission. During the prompt phase of LGRBs the spec- 
tral lag T^ is an ti-correlated with the peak luminosity a s 
shown by Fig. H (JNorris et al.l I2OOOI : lUkwatta etHlbOlOD: 
in co ntrast, short bursts have small t? ( Norris fc Bonnelll 
I2OO6I ) and occupy a separate area of the Lpk vs. r,^ /(I + z) 
parameter spac e (Fig. [6 1 and iGehrels et al.ll2006r ). Recently, 
iMargutti et al.l (120100 ) have demonstrated that, in strict 
analogy to their prompt pulses, LGRB X-ray flares deflne a 
lag- luminosity anti-correlation, where the lag is computed in 
the X-ray band (black dots of Fig. [6]). With the present work 
we complete the observational picture above, showing that 
flare candidates in SGRBs fall off the lag- luminosity relation 
defined by LGRBs: this furthermore supports a robust con- 
nection between prompt pulses and flare candidates in short 



^* Note that the same ~ 100 factor is found as the ratio of the 
isotropic energy emitted by LGRB flares and SGRB flare candi- 
dates: this is a direct consequence of their flare peak luminosity 
lO"'^, coupled to the very similar ui(tpk) 



ratio LSGRB/^lgrb 
evolution. 



"pk 
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bursts. At the same time this result points to some difTer- 
ences between LGRB flares and SGRB flare candidates. 

While SGRB p rompt p ulses arc significantly narrower 
than LGRB pulses l|Nakar et al. 2002 ; Nakar et al. 2002b :), 
fla res show instead com parable width in both classes (Fig. 
[2)). lHakkila et al.l (|2008r ) showed the existence of a lag- width 
correlation for prompt pulses: the wider the prompt pulse, 
the longer the lag. This b ehaviour has been recently ex- 
tended to LGRB flares bv iMargutti et al.1 (|2010b 'l. The sim- 
ilar width of LGRB flares and SGRB flare candidates im- 
plies that the lag-width relation cannot be invoked to explain 
the lag of SGRB flare candidates which are systematically 
shorter than expected from the lag-luminosity relation of 
LGRB flares. 

The physical cause of lags in the GRBs prompt 
emission is not yet understood: variations in the line-of- 
sight (ISal monson 2000); variations of the off-axis angle 



loka fc Nakamurall200 j) and rapid radiation cooling effects 
Schaefeill2004l ) are a few of the proposed models. In partic- 



ular it is at the moment unclear if lags in short bursts are 
small and non-measurable or intrinsically zero. According to 
the flrst possibility short and long GRBs would be powered 
by a similar, progenitor-independent physical mechanism, 
with SGRBs being faster evolving versions of LGRBs. The 
latter would instead point to some intrinsic differences. At 
present it is not possible to observationally distinguish be- 
tween the two scenarios r^l. 

In the case of flares, the situation is complicated by 
the limited and fixed energy band-passes used for the lag 
calculation (0.3-1 keV vs. 3-10 keV). The fundamental ori- 
gin of the observed lag is the spec tral evolution of a pulse 
(or fiare) profile to low er energies (jKocevski fc Liangj 120031 : 
iMargutti et al.ll2010bl ). The spectral peak energy Speak de- 
crease in time plays a major role in de termining the observed 
spectral evolution and lag value fsee iMargutti et al.ll2010bl 
for details). As a consequence, if the observed Speak does 
not cross the instrumental band pass during the emission, a 
limited spectral evolution will be detected and a short time 
lag determined. The shorter (when compared to what ex- 
pected from the LGRB flare lag-luminosity relation) time 
lag of SGRB fiare candidates might be partially a conse- 
quence of this observational effect (while possibly being in- 
trinsically larger) [fj: the results from Section 14.41 and 14.51 
suggest that long and short bursts basically differ in the in- 
trinsic time scale of central engine evolution (with SGRBs 
evolving faster). Since for LGRBs Sp °v f /-E!:°? p < 0-01 
(with Sp™F ~ 1 - 3 keV, observed value. IMargutti et al.1 
l2010bl ). the faster evolution of short bursts likely implies 
-EpSk!F/-EpSk!p < 0.01 for flares detected at the same ob- 
served tpk. This result translates into: Sp^^ p < 1 keV con- 



^^ From the prompt lag-width relation of iHakkila et al.l 1 1200811 : 
T^'LGRB/ LGRB ^ 0.0 1 - 0.1. Using w^°^^ ~ 0.05 s as typical 
value from lNakar et al.l l l2002bl ). we have t-t-S*^^^ ~ 0.5 - 5 ms 
assuming that the rf/wp ratio is universal. We typically resolve 
lags in SGRBs with a sensitivity of a few ms. This implies that 
we would be barely able to measure lags in the SGRB prompt 
pulses even if SGRBs were to follow the LGRB rT/top relation. 
^^ Note that the limited brightness of the flare candidates com- 
pared to the underlying X-ray continuum does not allow us to 
perform a one-to-one comparison with the spectral properties of 
LGRB flares to quantitatively check this possibility. 



sidering that Sp^S^^ ~ Sp^^ p as order of magnitude esti- 
mation (jGhirlanda et al.l 120091 ). According to this scenario, 
^pcak^ is below the XRT band for the majority of the emis- 
sion, possibly leading to a lag underestimation. The presence 
of this observational bias makes the interpretation of the en- 
tire lag-luminosity relation far from being straightforward. 
We stress that the dependence of the lag-luminosity on the 
choice of the flxed energy bands (both in the rest frame and 
in the observer frame) should be removed before addressing 
the physical interpretation of the anti-correlation. However, 
this topic goes beyond the scope of this paper and will be 
addressed in a future work. 

With this caveat in mind we note that if the energy 
Siso,F released by flares at different tpk is similar (as in- 
dicated by Section 14. 5p . then, considering that I/pk,F ~ 
Sjbo. f/w with the lag being positively correlated to the w 
(Ma rgutti et al.ll2010bl ). we would expect Lpk,F to be anti- 
correlated with t^^ as observed for LGRB flares of Fig. 
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In particular the normalisation 



iViag oc Siso.p, which implies A^iag"^ ~ lOOA^fl^^'' (since 
EIS^^ ~ 100£;?°f^, Table 0. This simple argument pre- 
dicts the SGRB flare candidates to be off the LGRB flare 
lag- luminosity relation of a factor > 100 as observed (the > 
inequality accounts for the possible underestimation of the 
real lag due to the fixed and limited energy band-passes bias 
of the previous paragraph) . This finding would support the 
presence of non-mesurable (but still non-zero) lags for the 
short burst prompt emission. 



4.7 The flare candidates internal origin 

The above indications point to a link between the properties 
of flare candidates and pr ompt pulses in SGRBs (for LGRBs 
this was demonstrated bv iMargutti et al.ll2010a) . This result 
would naturally arise if both kind of emission share a com- 
mon origin: we refer to this possibility as the internal origin 
interpretation. Theoretical models consistent with this pic- 
ture can be divided into two categories: models requiring a 
late-time GRB central engine activity and models where the 
central engine is not required to be active at the time of the 
flare detection. 

The second class of m odels includes the magnetic 
re-connection interpretation IjLvutikov fc BlandfordI l2003l : 
iGianniosI 12006) : flares would originate from residual, late- 
time magnetic re-connection events within the original out- 
flow (the same ejecta powered the prompt phase) triggered 
by its deceleration due to the sweeping-up of the external 
medium. The deceleration of the original ejecta during the 
afterglow phase causes an increase in the size of causally 
connected regions, thus enabling re-connection of increas- 
ingly larger magnetic structures. The smooth continuum 
would be instead afterglow emission from the shocked ex- 
ternal medium. 

Alternatively, flares and prompt pulses would automat- 
ically share a set of observational properties if they both 
directly arise from the GRB central engine activity (first 
class of models above) . If this is the case, the central engine 
would be active on much longer time-scales than previously 
thought (see e.g. lPerna et al.ll2006l : iRosswodlJOOTl : iLee et all 
:200S|); at the same time, the similarity of LGRB fiares and 
SGRB flare candidates in the I/pk/(L) vs. tpk/Tgo plane as 
well as in the AF/F vs. ipk/Tgo space would point to a 
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similar late time evolution of long and short GRB central 
engines. 

It is not possible to observationally discriminate be- 
tween the two scenarios using the present set of data. Careful 
modelling is required (Margutti et al. in preparation). 



5 CONCLUSIONS 

This work presents the first comprehensive analysis of flare 
candidates in SGRBs and provides a comparison to the ob- 
servational properties of X-ray flares in LGRBs with the aim 
to draw an observational picture of SGRB X-ray variability 
any theoretical model is required to explain. 

Our analysis shows that the origin of the SGRB X-ray 
light-curve variability is independent from the large-scale 
host galaxy environment and is not correlated to the de- 
tected afterglow life-time. We find that flare candidates ap- 
pear in different types of SGRB host galaxy environments 
and show no clear correlation with the X-ray afterglow life- 
time; flare candidates are detected both in SGRBs with a 
bright extended emission (EE) in the soft 7-rays and in 
SGRBs which do not show such component fSec. 14.2")) . We 
cannot exclude that flare candidates originate on top of faint 
(and undetected in the 7-rays) EEs. In particular, SGRB 
flare candidates are consistent with being drawn from the 
LGRB flare population when considering: 

1. The flare to prompt energy ration £'iso,F/£'iso,p (Fig. 
[51 Sec. 143)) : 

2. The relative variability time scale At/t < 1 (Fig. [T) 
Sec. [13}; 

3. The rest-frame flare width evolution with time i(;(tpk) 
(Fig.[l Sec. [13}. 

4. The hard-to-soft trend of the emitted radiation (see 
e.g. Fig.O. 

The main parameter driving the duration of the episodes 
of variability is the elapsed time from the explosion, with 
very limited dependence on the different progenitors, envi- 
ronments, life-times, prompt variability time scales and en- 
ergy budgets. The origin of the flare K;(tpk) relation must 
arise from what is in common for the long and short burst 
models. From another perspective this result implies that for 
tpk > 100 s the flare duration is likely to retain no memory 
of the variability time-scales of the original prompt mecha- 
nism at work. This would explain why the flare to prompt 
pulse width ratio is different for long and short GRBs. 

On the contrary, SGRB flare candidates significantly 
differ from the standard X-ray flare emission observed in 
LGRBs at the same tpk/(l + 2) in terms of: 

5. Peak luminosity L^gP^ _ o.OlL^g!^'' (Fig. [3] Fig. [1 



Sec. [431): 

6. Isotropic energy E^^^^^ ~ O.OIB;'^?^'' (Fig. [5) Sec. 

EM; 

7. Flare to prompt luminosity ratio Z/pk.p/ipk.p- Flare 
candidates in SGRBs are ~ 100 times dimmer than in 
LGRB; 

8. Relative variability flux AF/F (Fig: [7)): we find 
AF/F ~ 1 for all SGRB flare candidates (Sec. [43)): 

9. Lag-luminosity relation: like SGRB prompt pulses, 
flare candidates show shorter lags than expected from the 
lag-luminosity relation of LGRB flares (Fig. [6] Sec. l4.6p . 



However and more importantly, the differences listed at 
points 5., 6., 7. and 8. above disappear once the different 
time scale of evolution of the long and short GRB central 
engine as well as the different energy scaling of the two sys- 
tems are properly accounted for (Fig. [8] Sec. l4.5.l1 and [4753)) . 
This finding provides a connection between the properties of 
the detected SGRB X-ray light-curve variability and LGRB 
flares, suggesting a common, internal origin. As a result, 
we conclude that similar dissipation and/or emission mech- 
anisms are likely to be responsible for the prompt and flare 
emission in long and short GRBs, with SGRBs being less en- 
ergetic albeit faster evolving versions of the long category. 
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Figure B4. Same as Fig.[BT]for GRB 090607. 

APPENDIX A: SGRB FLARES: TABLE 
APPENDIX B: SGRB FLARES: PLOTS 



This paper has been typeset from a IJtjX/ X^TjrjX file prepared 
by the author. 
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Table Al. Best fitting parameters of SGRB flare candidates. From left to right: GRB, redshift; start time (ts), amplitude (A) and shape 
parameters (ti,T2) of the best fltting Norris et al. (2005) profile; width evaluated between 1/e intensity points; peak time (tpk); relative 
variability time-scale (Ai/i = ui/tpjj); relative variability flux AF/F: the value F is computed from the best flt of the continuous emission 
underlying the flare candidates (black solid line of Fig. [TJ; isotropic equivalent peak luminosity (Lpcak) a^nd energy (-Eiso) computed in 
the 0.3 — lOkeV observer frame energy band. 



GRB 
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Tl 
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AF/F 
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^pk.F 
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(s) 
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(s) 






(10-1^ erg s-1) 


(10«erg) 


050724 
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230.0 


12.29 ±0.88 
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0.26 ±0.03 
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Figure Bl. 0.3-10 keV count-rate light-curve of GRB 050724 and GRB 051210. Black solid line: continuous X-ray emission underlying 
the flare candidates computed as described in Section l2.2l dashed lines: best-fitting fiarc candidate emission; red solid line: best estimate 
of the total emission. The vertical dot-dashed lines mark the flare candidate onset times. Inset: Complete Swift-XRT light-curve. The 
yellow filled area marks the time window for the computation of the CCF lag. Middle panels: hardness ratio (HR) evolution with time; 
the HR is computed between 1.5-10 kcV (hard band) and 0.3-1.5 keV (soft band). Lower panels: Spectral photon index evolution with 
time as computed by Evans et al. (2010). 
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Figure B2. Same as Fig.|BT]for GRB 051227 and GRB 060313. 
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Figure B3. Same as Fig.|BT]for GRB070724A and GRB 071227. 



